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Abstract 

Motivated by its field theory interpretation, we study gravitational collapse of a min- 
imally coupled massless scalar field in Einstein gravity with a negative cosmological con- 
stant. After demonstrating the accuracy of the numerical algorithm for the questions we 
are interested in, we investigate various aspects of the apparent horizon formation. In 
particular, we study the time and radius of the apparent horizon formed as functions of 
the initial Gaussian profile for the scalar field. We comment on several aspects of the dual 
field theory picture. 



1 Introduction 



Understanding cooperative phenomena far from equilibrium poses one of the most challenging 
problems of present day physics. One source of methods is provided by the AdS/CFT corre- 
spondence [U El m H] which identifies a field theory without gravity with a string theory. In 
the low energy limit of string theory one can use the appropriate generalization of Einstein's 
equations to follow the evolution in time of the fields in the dual field theory. This is par- 
ticularly suitable for the description of strongly coupled field theories near a conformal fixed 
point. In this context the gauge/gravity duality opens a particularly important window: to 
study far from equilibrium phenomena one needs to study the evolution of Einstein's equation 
with appropriate boundary conditions. 

The AdS / CFT correspondence was readily expanded to include field theories at finite tem- 
perature [5]. In this context a field theory in equilibrium at finite temperature is dual to a 
black hole in asymptotically AdS spacetime. One very important development has been the 
establishment of the correspondence at the level of small fiuctuations. More precisely, the ap- 
plication of linear response theory for the near equilibrium region, id est, long wave lengths and 
low energies with local fiuid variables varying very slowly compared to microscopic scales (see 
[HI [TJ E] for reviews and complete list of references). The next frontier comes from the expec- 
tation that the evolution of spacetimes with the formation of black hole horizons is equivalent 
to non-equilibrium dynamics and evolution towards thermalization in field theory. 

Besides the general reasons to study far from equilibrium phenomena using the gauge/gravity 
correspondence one practical motivation comes from the RHIC and LHC experiments. Two 
crucial points are worth highlighting: (i) the observed Quark-Gluon Plasma (QGP) is strongly 
coupled and (ii) hydrodynamical description fits in a wide range of processes. Theoretical and 
experimental developments indicate that the QGP produced at RHIC is a strongly interacting 
liquid rather than the weakly interacting gas of quarks and gluons that was previously expected 
[H ini E] • A crucial piece of experimental evidence pointing to the need for methods applicable 
to strongly coupled theories comes from the fact that the produced plasma locally isotropizes 
over a time scale of Tiso < 1/m/c ( where c is the speed of light). The dynamics of such 
rapid isotropization in a far-from-equilibrium non-Abelian plasma can not be described with 
the standard methods of field theory or hydrodynamics [12] . In this manuscript we propose to 
study such rapid thermalization via its gravity dual - gravitational collapse. In this publication 
we largely expand, clarify and sharpen the results presented in our short letter |13j . 

String theory also provides a framework for understanding some superconductors holograph- 
ically [T3|IT5| (see reviews [ISIITT]). Moreover, some hard to understand properties of condensed 
matter system like the structure of non- Fermi liquids have been recently described using holo- 
graphic models [IB] [IS] (see [20] for a review). The question of thermalization is also crucial in 
condensed matter setups. There is an active interest in the understanding of the time evolution 
of a system following a quantum quench [211 122| [231 123]. A holographic approach to this area is 
showing promise in lower dimensions [25l [26] . The experimental observation of rapid thermal- 
ization in two-species mixtures of quantum degenerate Bose and Fermi Gases pTj constitutes 
another strong motivation for the study of thermalization using the gauge/gravity dictionary. 
The holographic description of these models always involves a black hole and the formation of 
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such a black hole will help our understanding of those systems at a more fundamental level. 

In this paper we present a full numerical analysis of the gravitational collapse of a massless 
minimally coupled scalar field in the presence of a negative cosmological constant in five space- 
time dimensions. We go beyond previous attempts within the AdS/CFT approach that relied 
on perturbation theory p8| or toy models for quench [291 SOI EI] . Our focus is on properties of 
the thermalization time. 

An intuitive explanation for the rapid thermalization time at RHIC within the context of 
the gauge/gravity correspondence is as follows. Local gauge- invariant operators are mostly 
sensitive to geometry near the boundary. As explained in [28], in the large limit which 
is the region explored by the gauge/gravity correspondence, trace factorization ensures that 
the expectation value of products equals the product of expectation values. Only one-point 
functions of gauge invariant operators survive. Other studies supporting rapid thermalization 
in the context of the gauge/gravity correspondence include [291 EH ESI 133] . 

The paper is organized as follows. In section [2] we set the stage for the problem, that is, 
we present the Ansatz and equations of motion. Section El contains details about the numerics 
involved and a brief discussion about the regime of very small amplitudes and gravitational 
collapse after a few bounces of the scalar field. Section HI contains our main new results from 
the gravity point of view where we present a systematic study of the time of formation of the 
apparent horizon. Since our work could be of interests to various communities (string theory, 
numerical relativity and non-equilibrium field theory), in section we survey the current 
literature for the benefit of the reader and explain how our results extend previous ones in the 
literature. We discuss various aspects of the field theory dual in section El We conclude in 
section El We also include a series of appendices where we collected a few facts that will help 
interpreting the numerical results presented in the main body of the paper; depending on the 
readers background they could be obvious. 



2 Collapse in asymptotically AdS^ spaces 

We consider the minimal Einstein-scalar field action with a negative cosmological constant 
A = 1)/(2L2): 
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The Einstein and Klein-Gordon equations of motion are respectively: 
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We tackle this system using numerical methods developed and tested in [3l], [35] which we will 
describe in more detail below. We choose the metric Ansatz to be 



ds^ 
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where a and a are each functions of only t and r. For convenience, let us define 



X = dr 



Y = -dt4>. 
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The Einstein equations can be conveniently written as: 
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Here we are using units where k = 1. The Klein-Gordon equation takes the form 



dtY = \dr (r^-X] -aa d^U. 
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The main difference of the system given by equations (12. 6[ 12.71 and 12. 9p with respect to the 
treatment presented in [S^ and [35] lies in the powers of r that appear in AdS^ there as opposed 
to AdS^ here. For example, in equation (12. 9p we have a cubic power of r rather than a square 
one. This higher power of r leads to a more singular behavior near r = which presents a 
numerical challenge. 

It is worth pointing out that equation (12. 8 p can be shown to be automatically satisfied when 
the other equations are satisfied. Numerically, that equation is used as a check to verify that 
the code is working properly. Namely the equation should be satisfied up to a numerical error 
due to the finite grid spacing. In the next section we will focus on this constraint in a fair 
amount of details. 



3 Comments on numerics 

We consider a massless scalar field, i.e., U = 0. We choose initial data of the form = 
y4exp(— (r — rof/a"^) where the amplitude A, center of the scalar field profile tq and profile 
width a are constants. This determines X through X = dr4> and we choose Y = X initially 
so that the wave starts out purely ingoing. We choose L = 1 and for the spatial grid to be 
< r < Tmax where r-^ax is a constant. For the simulations done in this work we use rmax = 10. 
Since in anti-deSitter spacetime, more resolution is needed at small r we do not choose the 
spatial points to be evenly spaced. Rather, for a simulation with N spatial points, the value of 
r at grid point z is r, = sinh(fc(z — 1)) where is a constant chosen so that = rmax- In the 
graphs below we use as a measure of the resolution of the simulation. 

For sufficiently large amplitudes the collapse process results in the formation of a black hole. 
This is signalled by the formation of a marginally outer trapped surface (also called an apparent 
horizon), where the outgoing null geodesies cease to diverge from each other and instead begin 
to converge. In spherically symmetric spacetimes, an apparent horizon has the property that 
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V^rVar = 0. In our coordinate system this would mean a^^ — )■ 0. The attentive reader whose 
intuition comes mostly from analytical solutions would expect that the horizon could as well 
be defined as gu which means a — near the horizon. That is not so in our situation 
as we stop the simulation based solely on the considerations of an apparent horizon discussed 
above and those considerations are independent of a. 

More precisely, our coordinate system breaks down when an apparent horizon forms, and 
so our signal that a black hole is forming is that a — )■ oo in the simulations. Actually, it suffices 
to set a fairly moderate value for the maximum allowed value of a, which we denote Omax and 
to stop the simulation whenever a reaches Omax, noting that a black hole has formed at that 
time. The position tah at which a = Omax gives us the size of the black hole and allows us to 
estimate its mass Mah = 3vr^(r^j|^ + r\^/L'^). Another notion of mass arises from the energy 
momentum tensor. This is the notion that best fits with the ADM mass, it is conserved, we 
denote it as asymptotic mass or total mass and can be computed as follows 





There is no a priori reason to expect Mah and Mbh to coincide as not all the matter density 
initially in the scalar participates in the formation of the apparent horizon. 

3.1 Understanding the numerical accuracy 

The code used in the simulations was originally used in a 4-d context in [3l], [35]. The precise 
version used here produced the results presented in [13]. 

In our previous work we tested that the code is working properly by considering the 
constraint quantity C = dta — raXY/3. Note that it follows from eqn. (12. 8p that C vanishes. 
However, the finite size of the grid spacing means that C will not vanish exactly in the simulation, 
but that instead in a properly working simulation, a smaller grid spacing should lead to a smaller 
C. In fig. (Il]) the result of a test of this sort is shown. We have plotted C for simulations with 
12800, 25600 and 51200 grid points. In all cases the simulations have A = 0.02, ro = 4.0 and 
a = 1.5 and the simulations are run to a time of t = 0.454. A gross estimate of the error 
directly from the figure suggest that error decreases by a factor of (1/2)^ pointing to a second 
order convergent code. 

In this paper we consider simulations with 12800, 25600 and 51200 grid points respectively. 
As a test of the accuracy we use a more invariant one. Namely, in fig. ([2]) we plot the L2 norm 
of the constraint as a function of time. The key feature that we want to point out is the value 
of the L2 norm of the constraint as a function of the grid resolution. The value itself is not 
the most important aspect but the way that it decreases as the number of points in the grid 
increases. 

There are a few observations that can be made directly from fig. ([2]). First, notice that 
we lose accuracy right at the time of formation of the apparent horizon corresponding to the 
almost vertical asymptote. What follows from the figure is that the simulation is very stable in 
determining the value of the apparent horizon as they all shoot up at the same value of time. A 
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Figure 1: The constraint quantity at various grid sizes. 
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Figure 3: Dependence of the asymptotic mass (see equation 14 .2^ as a function of the maximum 
value of a(r, t) which is used to determined apparent horizon formation. The asymptotic mass 
is expected to be constant and thus the plot can be understand as a test for the accuracy of 
the code. 



disastrous situation would have corresponded to losing accuracy at times that are significantly 
different from the apparent horizon formation time. It is also gratifying to see that the higher 
the resolution the closer we get to the time of the apparent horizon formation maintaining 
higher accuracy. We could, in principle, estimate the error in quoting values for times in our 
simulation based on the behavior of the L2 norms described in the figure. 

We have also performed more dynamical tests on the code. For example, we have checked 
stability with respect to the definition of the horizon amax by increasing its value tenfold and 
checking that the resulting radius and time of the formation of the apparent horizon are not 
affected to the precision we are working. In figure (|3]) we plot the value of the asymptotic 
mass for various values of the resolution as a function of amax- The precise definition of the 
asymptotic mass is given in equation f l4.2p : this quantity is constant in the simulation. What 
figure ([3]) shows is that for small values of amax the value of the mass is indeed constant as we 
increase amax we see that each simulation breaks but the point at which they break can postpone 
by increasing the resolution. In a more practical sense we are are following the maximum value 
of the function a{r,t), that is the plot can be also interpreted as related to the simulation 
time at which, due to the raising value of a{r, t) we lose accuracy. We also checked that the 
time and position of the apparent horizon formation were stable to the increment in resolution. 
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3.2 Oscillatory motion of the scalar field 

For collapse in asymptotically flat spacetimes, small amplitude initial data leads to a wave that 
is initially ingoing, then undergoes interference near the center and becomes an outgoing wave. 
For collapse in asymptotically AdS spacetimes, we would expect similar behavior, except that 
the outgoing wave should then reach the AdS boundary, bounce and become ingoing again, 
leading to another bounce, and so on. And indeed, this is what we find. Fig. shows the 
value of the scalar field at the center as a function of time. (0(0, t)). The parameters for this 
simulation are A = 0.0002, ro = 4.0 and a = 1.5. This simulation was done with 6400 grid 
points. Note that there are particular periods of time where the value of the scalar field at the 
center is non-negligible and that the middle of each such time period is separated from the next 
one by approximately vr. This is exactly what we would expect if the dynamics is mostly that of 
a massless scalar field on a background anti-deSitter spacetime. In the geometric optics limit, 
such a scalar field propagates along null geodesies, and for anti-deSitter spacetime (with L = 1) 
it takes a null geodesic a time of tt/2 to propagate from the center to infinity (see appendix [B] 
for details of geodesic calculations if needed.). 

0.25 
0.2 
0.15 
0.1 
0.05 
^ 
-0.05 
-0.1 
-0.15 
-0.2 
-0.25 

Figure 4: 0(0, t) for a small amplitude simulation. 

In our previous paper [13] we concluded that this oscillatory behavior could go on indefi- 
nitely. A more detailed study of this oscillatory behavior shows that the end outcome of this 
oscillatory behavior is black hole formation |36]. In a follow up note, Jalmuzna, Rostworowski 
and Bizon [27] evolved our initial data in a more accurate setup (allowing up to 2^'' + 1 grid 
points) and showed that a black hole forms after a time two orders of magnitude greater than 
those considered originally in [13]. An important proposal put forward in [36] states that a 
weak turbulent mechanism is responsible for the nonlinear instability of AdS with respect to 
arbitrarily small perturbations. The weak turbulent mechanism in question allows for a transfer 
of energy from larger to smaller scales. Corroboration for such mechanism was also recently 
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reported in |f38]. Notice that the incremental appearance of peaks in figure (jl]) is a visual 
counterpart to the mechanism proposed in 

In this manuscript we limit ourselves to show that we have explicitly verified that, in- 
deed, collapse is possible after a few bounces. We show, in fig. [5l data evolved for the 
A = 0.00084, a = 1.5, ro = 4 with a resolution of = 12800 points. We used [37] to find 
the most appropriate range of parameters for our purpose. Note that the scalar field bounces 
once and then gravitationally collapses, we show that in the figure by tracking the value of 
the metric function a{r,t) near the position of eventual horizon formation. In the simulation a 
black hole forms at vah = 0.0049 after a simulation time of tsiM = 4.568. 




SIM 



Figure 5: a{r ^ r^H^t) for a small amplitude simulation that bounces once and then gravita- 
tionally collapses. 



4 Black hole formation 

The quantity that we are most interested in is how long the process of black hole formation 
takes, i. e., we would like to know how much time elapses between the initial time and the time 
a marginally outer trapped surface forms. For simulations it is natural to use the coordinate 
t as a measure of time and we have used it previously in [13]. Here we continue to use that 
simplified notion of time, a more precise definition of the field theory time is described in 
appendix ([C]) and amounts to correcting the coordinate time to incorporate the right field 
theoretic normalization. We denote by ^ah the simulation time t at which an apparent horizon 
first forms. Note that our coordinate time (and similarly the field theory time) t has geometric 
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meaning in that it is the coordinate that is orthogonal to the area coordinate r. We discuss 
various relevant notions of time in more detail in appendix (lUj) . 

We also want to know how the time of apparent horizon formation depends on the initial 
amplitude of the pulse A. This dependence is shown in fig. which gives ^ah for several 
simulations with different values of A. For each of these simulations we have a = 1.5 and 
To = 4.0. These simulations were done with 51200 grid points. Similarly, in figure ([7]) we 
plotted the dependence of the radius of the apparent horizon as a function of the amplitude. 
The masses of the black holes formed in the simulations range from 23.5 for the smallest 
amplitude to 8.52 x 10^ for the largest amplitude, that is, about three orders of magnitude. 




Amplitude Amplitude 



Figure 6: The dependence of ^ah on A Figure 7: The dependence of tah on A 

Some aspects of figures ([6]) and ([7]) can be understood using the properties of null geodesies 
in anti-deSitter spacetime. We are going to slowly build up our intuition on this account. 

First, fig. d?]) shows that the radius of the apparent horizon grows with the amplitude A. 
This is intuitively clear as the higher the amplitude A the higher the initial energy stored in 
the scalar field and consequently the higher the value of the asymptotic mass Mbh defined in 
equation ( 13. ip . 

To gain a better understanding of the time of black hole formation tAH we also study its 
dependence on the initial width of the pulse a. This dependence is shown in fig. (jH]) which gives 
^AH for about a hundred simulations with different values of a. For each of these simulations 
we have A = 0.02 and tq = 4.0. These simulations were done with a radial resolution of 51200 
points. 

The main intuition that follows from this graphs is that once the width has increased to a 
given value, the time of formation of the apparent horizon remains approximately constant. It 
is instructive to understand this result in a semi-analytic fashion. Rewriting the metric function 

a as 

, r2 M(r,t) , , 
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we have that equation (I2.6p takes the following form for the local mass M(r, t) 

d 



dr 



M{r,t) 



IT 
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(4.2) 



The asymptotic mass called Mtotai in equation (13.11) is simply the limit as r — > Vmax of M(r, t). 
We can estimate, in the thin-shell approximation, the dependence of M on the width of the 
initial Gaussian profile. Namely, in the thin-shell approximation (cr ^ tq) we evaluate the 
above integral at t = where the X = Y and we can approximate a ~ 1/r, then the leading 
term in a is 

o(^] (4.3) 



M 



^' 2 a \ro. 

This behavior of the mass leads to a dependence of the form r^n ~ a~^^^ which is precisely 
what we see in figure for small a, in the range a > 1 other contributions to M flatten the 
curve. 




width 




Figure 8: The dependence of ^ah on a 



Figure 9: The dependence of tah on a. 



This figure allows us to conclude that for a given amplitude A increasing the width a in the 
original pulse amounts to decreasing the mass of the black hole that might eventually form. 

The time it takes a null geodesic to travel from the center of the Gaussian profile at Tq 
to the apparent horizon radius tah can be explicitly calculated in AdS and the result is (see 
appendix [Bl) 



tAH 



^arctan (^-j- 



arctan (^^~" j 



(4.4) 



One of the main results of our paper is contained in the figure f|TO|) . The figure contains the 
result of about a hundred simulations, we have plotted t^H versus r^H from the simulations. In 
the same figure we have superimpose the analytic expression for the time it takes a null geodesic 
to go from the center of the original Gaussian packet to the horizon radius. We have also 
presented a best fit with an analytic form that closely follows the result of the geometric 
approximation. 
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Figure 10: The dependence of ^ah on the radius of the apparent horizon 
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The intuition that the time of horizon formation would be closely determined by the null 
geodesic time was a prevalent intuition in the literature. Our simulations show that the time of 
formation of an apparent horizon is still of the same order but slightly bigger than the one given 
by a null geodesic as figure (fTOj) nicely shows. Clearly, the width in the original Gaussian is 
playing a prominent role. Since we are solving the full nonlinear system of Einstein and Klein- 
Gordon equations we are able to take into consideration configurations far from the thin-shell 
limit. The full understanding of the shape of the Gaussian profile in the time of formation of 
the apparent horizon will require more work. We will focus on the prevalent role of the width 
in what follows. 

We finish this section by advancing an interpretation of the best fit to the data as presented 
in figure ( fTOj) . First, the value 1.3751 turns out to be well approximated by arctan(5), this 
suggests that the effective center of the Gaussian profile closer this value than to ro = 4. This 
is a direct consequence of the properties of AdS. Second, we interpret the coefficient of 1.0441 
as the extra time that it takes a shell of finite width to collapse. Thus figure ( ITO|) with its 
best fit curve can be understood as summarizing the effects of a null geodesic falling from an 
effective distance of ro = 5 in AdS and also the effect of the extra time that takes the Gaussian 
to shrink. 



4.1 A quantitative picture of gravitational collapse 

With the hope of clarifying the main properties of the space-time during the gravitational 
collapse we present fig. flTT]) . In the figure we plot the Klein-Gordon energy density. Consider 
a unit timelike vector n = a^^d/dt, then the energy density is defined as (see appendix R|) 

p = T^.n^n" 

f 6 »2 

Based on this energy density we will plot a slightly different quantity which focuses on the 
scalar energy 

PNormir, t) = ^ (4.6) 

/ ^(x^+n 



The advantage of this quantity is that for each value of time we normalize the energy to its 
total value on that line. This is simply a choice of presentation, otherwise, it becomes hard to 
follow what regions are really dense at a given time. 

The simulation presented in fig. (11 II) was conducted for A = 0.125, o" = 1.5, ro = 4 and with 
N = 51200 grid points. In fig. ( ITTj) the vertical axis represents time and the horizontal axis 
represents position r. There are various properties to note about our measure of energy density 
for the scalar field. The energy density is proportional to the derivative of the scalar field 
and, from the initial condition, derivative in time and in the radial directions are correlated. 
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therefore, we see a dark blue spot at r = tq = 4 corresponding to a zero of the energy density 
at t = 0. One might expect a symmetry around that point, however the expression for energy 
contains an effective extra power of and therefore shifts the maximum value to r = 5. The 
black vertical line in figure ( ITT]) represents the spatial position where the apparent horizon 
eventually forms. The horizon only actually forms at a later time in the simulation but we 
draw this spatial black vertical line to highlight the position. To guide the reader we have also 
drawn two white dashed lines between which 90% of the Klein-Gordon energy density of the 
scalar field is contained. 

What we see clearly is that, to a good approximation, the front part of the scalar wave 
crosses the apparent horizon and the second part moves effectively in a Schwarzschild-AdS-like 
background and experiences a marked stretching. 
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Figure 11: Time evolution and radial distribution of the energy density stored in the scalar 
field with Gaussian profile as the apparent horizon forms. 

Let us also remark in this section a fact that will be important in the field theory interpreta- 
tion of our results. The deep blue regions in figure (ITT]) are regions where the energy density of 
the scalar field is basically zero. The equations of motion, or more generally Birkhoff 's theorem, 
imply that the spacetimes in those regions can only be AdS^ or the Schwarzschild black hole in 
AdS^. Therefore a boundary configuration that explores only the asymptotic region will only 
detect AdS^ or the Schwarzschild black hole in AdS^. 

For completeness, in fig. f|T2|) . we also present the energy density graph for a situation where 
the scalar field bounces once and then collapses. The parameters of the simulation are exactly 
those presented in the previous section. Namely, A = 0.00084, a = 1.5, ro = 4 with a resolution 
of = 12800 points. It is worth noticing that we have performed the same calculation with 
other resolutions and the results are indistinguishable from those presented here. 
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Figure 12: Time evolution and radial distribution of the energy density stored in the scalar 
field with Gaussian profile as the apparent horizon forms at a very small radius tah = 0.0049 
after a bounce of the scalar field. 

4.2 Comparison with previous results in the Uterature 

The question of gravitational collapse in asymptotically AdS spaces has attracted a lot of 
attention largely due to its application in the AdS/CFT correspondence. Some of the first 
works, for example [391 SOI SI] discussed the setup and its implication in the context of the 
gauge/gravity correspondence. An interesting model for collapsing matter shells in analogy with 
the Oppenheimer-Snyder solution was discussed in [42]. Other versions of the Oppenheimer- 
Snyder model in asymptotically AdS spaces were also considered in [i3l HI] . In the numerical 
relativity community, similar problems have received a lot of attention. For example, collapse 
in asymptotically AdS^ has been addressed with beautiful results in |l5l HHl HTJ HH] . Similarly, 
some work in AdS4^ was presented in, for example, [19]. More general dimensions have been 
briefiy considered in [50j. Aspects of potential cosmic censorship violations in asymptotically 
AdS4^ collapse have also been addressed by Garfinkle [HI EH]. Recently, a discussion of black 
holes in a box, which is an important step in the full construction due to the fact that AdS is not 
globally hyperbolic, was presented in [511 [52] . More recently the question of turbulent instability 
in asymptotically AdS spaces has been studied in the spherically symmetric [361 EZ] cind non- 
spherically symmetric [38] setups; the complete set of implications for the AdS/CFT remains a 
truly interesting problem. There are other approaches that have been pursued with similar field 
theoretic intent. For example, a different way of setting up the problem is due to Chesler and 
Yaffe [531 El] who considered a time dependent boundary metric (7^^ which creates gravitational 
radiation that propagates from the boundary to the bulk. This gravitational radiation leads to 
the formation of a horizon. The anisotropic nature of the problem is particularly relevant for 



14 



applications to collision situations but also for the quantum quench setup discussed recently 
directly in field theory [21], [22], [55], [5S], [23], [21] ■ Another example of a different approach 
has been presented in a beautiful recent paper [57] where the authors studied the approach to 
the transition to the hydrodynamic regime using newly developed techniques for higher order 
hydrodynamics in the context of the boost-invariant plasma originally proposed in [58] . Namely, 
the authors of [5Z] studied the far from equilibrium dynamics of a system with boost-invariant 
symmetry, arguably the dual of Bjorken fiow [59j in field theory; they showed that using higher 
order hydrodynamics the approach to the linear hydrodynamic regime occurs in time scales 
inversely proportional to the equilibrium temperature. 

The thin Vaidya shell approximation - Of the many setups discussed in the htera- 
ture, the one that can be most directly compared to our work here is the one presented in [28] . 
This work provided, to our knowledge, the first systematic analysis of the thin shell Vaidya 
collapse as the dual of rapid thermalization in field theory. 

Let us briefiy discuss the setup in [28] which considered an action of the form (12.11) . Their 
metric is spherically symmetric and is given in Vaidya type coordinates 

ds'^ = 2drdv - g{r,v)dv'^ + f{r,v)dn'^ 
(j) = (f){r,v). (4.7) 

where dQ'^ is the line element of the unit S^. The crucial physical information is stored in 
the scalar field profile which is (poi^) < e for < f < T and otherwise vanishes. In [21], a 
perturbation theory in e was developed. The situation we treat is not precisely the same as 
that of [2H]: our initial data is different, and the fact that we use different coordinate systems 
makes direct comparisons somewhat involved. Nonetheless, there are many similarities. Our 
Gaussian profile falls off so fast that it might as well be of compact support, and our choice 
of ingoing waves means that at large r the profile is essentially a function of an ingoing null 
coordinate v. We can also choose our parameters A and a (corresponding respectively to the e 
and T of [21]) sufficiently small to be within the perturbative regime of [21] or sufficiently large 
that that perturbative regime is no longer valid. Furthermore, in both cases many of the main 
features of the collapse process seem to depend simply on the approximate propagation of the 
scalar field along null geodesies until the shell becomes sufficiently small that a trapped surface 
can form. The main advantage of our method is that we can state rather precisely such things 
as when a black hole forms and how large it is without having to be in the regime where the 
perturbation expansion of [2H] is well approximated by its first couple of terms. 

To make contact with the Vaidya form of the metric we consider, in the metric 12. 4[ the 
following coordinate transformation: 

t = fi{r, v), r = r. (4.8) 
The condition that the resulting metric has v = const, as ingoing null geodesies leads to 

?t = (4.9) 

or a 

We can adjust our metric explicitly to contain a null incoming coordinate: 

ds'^ = —aldv'^ + 2a^a„(iw dr + r'^dil^, (4-10) 
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where = [dn/dv) a and = a{r,v). One can now specialize the general form above to the 
particular case of the Vaidya metric in global coordinates which is 

ds^ = -(1 + _ ^^)dv^ + 2drdv + r^dnl (4.11) 

One important property of this metric which justifies is wide range of applications is that it 
solves Einstein's equation with null dust. Dust is characterized only by its energy density which 
needs to be 

In various recent works in the literature |25l |26l [291 EI] M{y) ~ tanh(i;). It is possible to 
look at figure ( ITTi) and see how the collapse picture differs from the thin shell model described by 
the Vaidya metric. It is worth noticing that figure ( ITTi) is rather representative of the hundred 
or so values of the width a that we have considered in this paper; only the particular details 
change but the qualitative picture is fairly universal. 



5 Comments on the dual field theory thermalization 
5.1 Field theory setup 

Our field theory setup is rather specific and we clarify it now. The overall field theoretic 
question one is interested in answering is: How does a field theory react to a rapid injection 
of energy? This is precisely what the RHIC and LHC (ALICE) experiment are all about for 
large values of the energy of the colliding particles - How does QCD matter behave under such 
collision? Now the collision is clearly anisotropic as one has two gold atoms colliding. Direct 
head-on collisions can potentially allow for semi isotropic treatments. A gravity approximation 
to the RHIC and ALICE experiment has been developed in the context of numerical relativity 
by Chesler and Yaffe in |53l |54l |60]. Several authors have considered simpler versions of this 
problem by utilizing colliding shock waves but we do not discuss those works here. A natural 
time scale in such collision process is the isotropization time which in the case of RHIC is 
Tiso < Ifm/c [12]. 

It is worth pointing out that in the presence of an isotropization time scale we have to come 
to terms with the existence of another scale which is the time scale where the hydrodynamic 
regime becomes applicable. Recall that the thermalization time is, by definition, the time after 
which particles follow a Boltzmann distribution. The hydrodynamic time, that is the time 
when the field theory is well described by the hydrodynamic approximation, can be applicable 
either after the thermalization time or when the pressures in the rest frame become equal 
(isotropization). In the interesting paper of the Chesler- Yaffe series |60] evidence is found for 
the hydrodynamic regime to be applicable even when the difference between the pressures in 
the longitudinal and transversal direction is still large, that is, the hydrodynamic time scale is 
smaller than the isotropization tim^. Note that other approaches, for example [57], tend to 
identify the hydrodynamic and isotropization times. 

^We thank D. Mateos for this important clarification. 
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Our setup has one major difference. Namely, our injection of energy is spatially homoge- 
neous, that is, it is the same in all points of the field theory space at a given time. We thus study 
a situation that is similar but not exactly equal to a quench. In this sense our setup is perhaps 
less directly applicable to RHIC-type experiments of collision but rather speaks directly about 
more universal properties of strongly coupled field theories at large which are dual to gravity 
theories. What we set out to study in this paper, using gravity methods, is the time in field 
theory between the injection of energy at t = and the formation of the quark gluon plasma. 
We call this time the thermalization time. Of course, other interesting time scales might be 
present in a given experimental setup. One example is the above mentioned isotropization time, 
another natural time is the time at which the hydrodynamical approximation becomes a valid 
description of the quark gluon plasma. The thermalization time, as defined in our context, is 
expected to be of the same order but slightly smaller than these other two. 

To conclude this part by emphasizing that there are various time scales in the problem some 
of which we list: 

• Hydrodynamic time scale: scale after which the liquid is well described by hydrodynamics, 
small velocities, long wavelengths. 

• Isotropization time: The time scale after which the pressure differences in various direc- 
tions is small. 

• Thermalization time: Time after which the particle distribution follows the Boltzmann 
distribution. 

5.2 Field theory regimes 

As we know, in field theory not every collision or injection of energy will lead to formation of 
a thermal state. Experiments that involve low energies, and low number of particles will be 
better described as scattering states. As the energy and number of particles involved increases 
we expect the final state to be better approximated by a thermal state. As demonstrated in 
the RHIC experiment, at strong coupling the energy needed to generate quark-gluon plasma is 
lower than the naively expected, however still finite. 

Recent results due to Bizon and collaborators [361 EZ] and also corroborated in [38] indicate 
that any perturbation in AdS ultimately leads to gravitational collapse. We would like to 
explain how these results are reconciled with the standard field theory expectations. 

From the gravity point of view Bizon's result translates into a formation of a thermal state 
for a collision of arbitrary small energies or for very small amount of energy. However, when 
translated into a field theory language we need to incorporate the large Nc scaling. The energy 
scale for the collapsing shell involve Newton's constant k^, when properly translated into field 
theory language we counter a N"^ scaling of all energies, equivalently an order of A*"^ particles. 
In terms of the AdS/CFT duality the gravity approximation is valid only for infinite Nc, thus 
in term of field theory we always have enough particles to form a proper thermal state. If we 
were to discuss field theory collisions that are too small to generate quark-qluon plasma, we 
should introduce stringy objects and quantum excitations in the AdS side. 
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One can summarize the above discussion by identifying two regimes: 



• Regime without backreaction: The field theory description of a colhsion with very low 
energy or a low energy quench should find its gravity description in the geodesic motion 
of objects in AdS. 

• Regime with backreaction: The formation of a black hole with nonzero area in the gravity 
side is equivalent to the formation of a quark gluon plasma in the field theory. The fact 
that this always happens when considering a scalar field in Einstein gravity with a negative 
cosmological constant is a consequence of the large Nc limit. We will show more rigorously 
in the next section that this scaling does indeed take place. 

5.3 Hawking-Page transition and gravitational collapse: Canonical 
and microcanonical ensembles 

There are many arguments leading to the fact that the gravity calculations presented here 
pertain to the microcanonical ensemble. This questions has been revisited several times in the 
context of general relativity. First, conservation of energy tells us that we are in a system with 
a fixed amount of total energy. Further, since black holes in fiat space have negative specific 
heat they can cool themselves by transferring energy to the surrounding and consequently be 
able to transfer more energy. We expect that small black holes in AdS will behave similarly. 

From the field theory point of view sometimes the canonical ensemble is relevant in which 
case it is worth revisiting the role of the Hawking-Page phase transition in the context of the 
AdS/CFT correspondence. The interpretation of the Hawking-Page transition in the context 
of the gauge gravity correspondence leads to the result that below a certain temperature, that 
is, below a certain radius of the black hole, the background that is relevant in the thermody- 
namic competition is no longer the Schwarzschild black hole but rather thermal AdS [5j . This 
result means that for low temperatures the field theory goes into a confined phase rather than 
remaining in the quark gluon plasma. These suggests that the very low temperature black holes 
are not relevant from the field theory point of view as they correspond to an unstable phase. 

The Hawking-Page transition occurs, using Witten's notation [5J, at r+ = L. Namely 

_ vol(g"-^) L'rl-' - rl+^ 

4Gn nrl + {n-2)L^ ^ ^ 

For r_|_ < L the thermodynamical competition is won by thermal AdSn+i] for r_|_ > L the 
dominant themodynamical phase is the black hole. In the work [13] and the discussions of 
[36| [37] . black holes with r+ < L = 1 should not be considered as part of the description 
of a dual field theory since in that range the field theory stable phase is best described by 
thermal AdS. A counting of the degrees of freedoms allows to identify this transition with the 
confinement /deconfinement transition given that the entropy of the black hole is of the order 
while the entropy of thermal AdS is of the order iV°. 
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5.4 Boundary stress energy tensor 



Translating gravity quantities into field theory terms is, in general, a task that has to be dealt 
with on a case by case basis. In this section we will simply use the dictionary for the canonical 
case of the duahty in the context of classical strings in AdS^ x dual to A/" = 4 SYM: 

^ L^ = 47ig,a'\ (5.2) 



k2 647r4c/2«'4' 

For many quantities, for example Kubo-type formulas, one can verify that the gravity quantities 
combine in a way that there is no dependence in a' leading to 

Our goal in this section is to extract the field theory stress energy tensor. The prescription for 
extracting the boundary stress energy tensor using bulk gravity data is well established in the 
context of the AdS/CFT correspondence [SI][S2] and reads: 



(T,,) = -^\imr%,, (5.4) 

47r r— >-oo 

where fi, v run over the field theory directions and ^g^v is defined as the difference of the resulting 
metric from the boundary metric corresponding to AdS. We are interested in the boundary 
or field theory energy which we relate to Too in the sense above. We are thus, interested 
in studying the asymptotic behavior of the function a(t,r) from which we read bg^^ above. 
The most efficient approach is to approximate, in the asymptotic region, the function a by 
its Schwarzschild-AdS expression. Using the equations of motion we could, equally effectively 
study how well approximated is the function a by its Schwarzschild-AdS expression which is 
what we actually do in the simulations. 

Figure [13] shows how well the function a{r,t) is approximated by the Schwarzschild-AdS 
expression at times close to the time of formation of an apparent horizon. To aid understanding 
the fit better we have plotted the difference between the function a(r, t) and its Schwarzschild- 
AdS expression. It is remarkable that this particular values of the simulation parameters the 
difference is better than 10~® as we approach the asymptotic regions. 

Finally we consider the same plot for a hundred different values of the amplitude of the scalar 
field. The fitting is better than half a percent for the majority of the values tested. Except 
for the two anomalous points, larger amplitudes yield a better fitting to the Schwarzschild-AdS 
form of the metric. Notice that at small horizon radius we do expect a deviation between Tbh 
and r^H as at the time of formation of the trapped surface, some portion of the scalar energy 
is still un-trapped. 



5.5 RHIC and LHC isotropization times 

In this section, despite all the disclaimers made previously, we take our results in the direction of 
application to actual experiments. Let us remark again that we understand that thermalization 
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Figure 13: The difference close to ^ah of the metric function a{r,t ~ tAn) directly from the 
simulation and a fit with the static Schwarzschild-AdS metric. 



Relative error between r„,, and the AdS-bh horizon radius r^^ 

AH bh 




Figure 14: Relative error in determining the horizon radius from the simulation value of the 
apparent horizon and from the fit of the metric functions according to the static metric of 
Schwarzschild-AdS 
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time and tAu could be widely different. Our goal is not to make a concrete numerical prediction, 
rather, we aim at simply illustrating where our results might fit. One of the results of the gravity 
simulation is a description of the dependence of the time of formation of an apparent horizon as 
a function of the black hole mass. When translated into field theory terms we could argue that 
this graph approximates the dependence of the thermalization time (formation of the quark 
gluon plasma roughly) as a function of the energy. Below, in figure [13 we plot such graph. 

0.9 
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0.3 
0.2 
0.1 



Figure 15: The dependence of ^ah on the mass of the formed black hole 

The best fit to the curve in figure f|T5l) could be understood as a curve describing tThermaHzation{E) . 
We shall loosely interpret E as the energy injected into the field theory. We consider explicitly 
a fit for large masses which clearly follows from the geodesic approximation presented in fig. 
( |T5l) . Note that this expression is a deviation of the naive one that can be obtained in the null 
geodesic approximation for AdS. Namely, from equation (14. 4p it follows, in the large mass limit 
there, that t ^ _£;-o-25 where we relate E to the radius of the black hole using the Schwarzschild 
solution. The precise numbers appearing in fig. (fTSll have the same interpretation discussed in 
section HI Without providing further evidence we mention that we have also considered similar 
graphs for different widths of the scalar field and the results come out in the same range of 
the exponential. We do not claim that this is a universal behavior but within our modest 
context this exponent has come out about the same for different values of the widths. It will 
be interesting to explore this potential universality in the Iah versus Mbh graph further. 

Hoping that the analysis we performed is universal enough we dare to compare our results 
with experimental results. There is really no reason this might work as we are comparing times 
in very different experiments. As explained in the introduction the thermalization time we 
compute here and the isotropization time at RHIC and LHC correspond to different physical 
situations. However, there is the possibility that the time we compute here is of the same order 
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as the isotropization time as it describes a fairly universal property of the dual field theory, 
namely, its response to a a sudden influx of energy. With that disclaimer, we could conclude 
that ^^^"^ = '^tfj^""^} . We use [63] to estimate that the central Au-Au collisions occur at 
approximately y^s^viv = 0.2TeV while the data of the ALICE collaboration at CERN [M], [65] 
suggest that the central Pb-Pb collisions at ^/sNN = 2.76TeV. Plugging in the numbers we get 
tlhc ~ -^trhic- 



6 Conclusions 

In this paper we have considered gravitational collapse of a scalar field coupled to Einstein 
gravity with a negative cosmo logical constant. We have focused on the dependence of the 
time of formation of an apparent horizon as a function of the parameters characterizing an 
initial Gaussian profile for the scalar field. We varied the width and amplitude of the initial 
Gaussian profile and studied how it affects the time of formation of an apparent horizon. We 
obtained results for masses ranging over two orders of magnitude. In a more detailed study we 
verify that the dependence on the width is milder than the dependence on the amplitude. Our 
simulations indicate a saturation of sorts, in the time of formation of an apparent horizon, as 
we increase the width of the Gaussian profile. In conjunction with [13j we continue to build 
on the foundation that will allow us to provide, in full detail, various other properties of the 
thermalization process by studying its gravitational dual. We have proposed that these studies 
might provide rough estimates for quantities that depend on the thermalization time. Even 
though the concept of thermalization itself depends on the precise set up, gravity studies should 
provide a rough estimate. 

We conclude enumerating some of the directions that we believe would be interesting to 
pursue. 

• We have verified that for some parameters in the simulations the scalar field does collapse 
after a few bounces as explained in [35] and ^7\. The relevance of this result for the 
microcanonical ensemble in the field theory dual is hard to overstate and should be 
clarified. 

• What is the role of the potential in the scalar field? It is our expectation that considering 
a potential term of the form m^0^ could significantly impact our results as an underlying 
intuition throughout the paper has been the behavior of null geodesies. 

• Spherically symmetric approximation is poor for describing collisions. Key properties, 
such as elliptic flow, are missing. Although the techniques involved in going beyond the 
spherically symmetric framework are substantial it is worth pursuing. 

• Other interesting time-dependent phenomena have recently been tackled using gauge/gravity 
correspondence. In particular aging phenomena have been discussed in [BSl EZl EB]- It 
would be nice to have a more precise gravity description along the lines obtained in the 
context of this time-dependent far-from-equilibrium phenomena, aging. Note that [68] 
obtained the field theory correlators using a gravity background and the corresponding 
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gauge/gravity prescription for correlators. One really intriguing question would be to 
understand precisely what initial quenches lead to black hole formation and which lead 

to an aging situation. A numerical realization along the lines presented in this paper of 
aging phenomena would be a great insight into the structure of aging phenomena. 

• Finally and perhaps most interesting conceptually: Can we use our solutions to set up 
a laboratory to study the black hole information paradox? There is no doubt that the 
information paradox is resolved in the context of the gauge/gravity correspondence as we 
have a map between a gravity theory with black hole formation and evaporation and a 
unitary field theory. The answer to the information paradox has rather become a question 
of how it gets resolved in our models of gravitational collapse. A detailed study of black 
hole formation with the various stringy probes and correlation functions is a tantalizing 
route to understanding this question better and to an eventual explicit resolution of some 
of the most interesting black hole puzzles. 
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A Equations of Motion 

The Einstein tensor Gf^i, — R^^ — ^Q/^uR is 



2 r 



G 



3a 

ra 
3 

a 

C - 1 



a(l — a^) + rc^rO 



'1 — cc^)a + r^rO; 



— ra^dra{2a + S^a) + aa^{a + r(29j.Q; + rd^a) + r^a^{dtadta — ad^a) 



Gee = Gyy sin^ x, 



G^^ sin^ X^i^^ ^- 



(A.l) 



The stress energy tensor is written in a form that contains the contribution from the cosmo- 
logical constant 



T 
Ttt 

Tfr 

T 



XX 



dt(f)dr(f), 



(A.2) 



B Geodesies in global AdS^ 

In this appendix we recall a few known facts about AdS in global coordinates. We discuss a 
massless scalar field that should move along a null geodesic. Recall that a null geodesic reaches 
spatial infinity in finite asymptotic time. Namely, considering the metric 



ds' = -\\^ — \dt^ + ^ + r^dVil. (B.l) 

12 



For a radial null geodic we have 

d\ ' d\ 1 + 



* = ±E, * = (B.2) 



L2 
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In particular such a null geodesic reaches r — ?► oo starting at r = in 



II. (B.3) 



In our analysis of gravitational collapse we will naturally observe that for solutions where the 
amplitude of the scalar field is not sufficient to form a black hole the wave will roughly oscillate 
many times which is a combination of reaching asymptotic infinity in finite coordinate time 
and the fact that we impose Dirichlet boundary conditions at infinity since AdS is not globally 
hyperbolic. 

The situation for massive geodesies is similar. Starting from the metric (IB.ip . the relevant 
equations of motion are 

dt E ^ r"^ r, , 



The solution of this massive geodesic is 



1 



r(A) = LV^2 _ isin _ (B.5) 

Interestingly, the time that it takes this particle to reach its maximum radius Vmax = L\/ — 1 
is precisely as in the case of the null particle 

E 71 

dr =-L. (B.6) 



C Times: simulation, coordinate, physical and field the- 
ory times 

In this appendix we clarify the different measures of time that we used throughout the paper. 
The coordinate time is simply t. We make sure that at each step the simulation time step dt is 
at least smaller than the corresponding physical time step associate with an increment in the 
radial direction by dr, that is, dt < {a/ a) dr. In the bulk of the paper we use this coordinate 
or simulation time and call it tsim- 

The metric given in equation (12. 4 p has still a freedom that we would like to fix. Namely, the 
form of the metric is unchanged if we apply a transformation of the form t — )■ /(t), a — )■ af'{t). 
We would like to use a gauge in which the time at infinity is normalized in the standard form. 
However, that was not the choice we took in the simulations, rather we used a{r = 0,t) = 1. 
Because of this choice we use in some of the graphs what we call field theory time 

dtpT = a{r = r^nax, tsim) a{r = r^ax, tsim)dtsim- (C.l) 

In the beginning of the simulation the difference between tsim and tpT is essentially a constant 
that accounts for the right field theory normalization. However, near the time of formation of 
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Figure 16: Plot of the field theory time tpT versus the coordinate or simulation time tsi: 



the apparent horizon the situation is quite different as can be seen from the plots in figure ( !T6|) . 
Note that that figure ( fT6l) provides a check of the stability of our definition as a function of the 
grid resolution. Namely, despite the violent dependence around the formation of the apparent 
horizon different resolutions behave remarkably close to each other are are essentially indis- 
tinguishable in the plots. Since the relation between the field theory time and the simulation 
time is via an integration, figure ( !T6|) shows that the accumulated error is negligible among the 
various grid resolutions that we used. However, near the time of formation of the apparent 
horizon tpT breaks down and should not be used to define field theory thermalization time. 

To develop a better intuition of the implications of using field theory time, we consider the 
constraint plotted as a function of the field theory time. Notice that the constraints have a 
milder behavior near the time of horizon formation as can be seen in the figure f|T7|) . This figure 
should be compared directly with figure ([2]) of section 13.11 The key feature is the smoothing 
around the time of formation of the apparent horizon which indicates that the simulations are 
not really losing accuracy as one could infer from figure ([2]) of section 13.11 

Motivated by the remarks of section (14. 2 p we attempted a computation of the ingoing null 
geodesic time denoted by v. This coordinate is essentially what would be required to go to the 
Eddington-Finkelstein coordinates. This null time f is a key ingredient in various AdS/CFT 
setups, certainly in computations of spectral functions it is this time that naturally fixes the 
boundary conditions near the horizon. It also figures prominently in the hydrodynamic context 
(see O O E]). The setup used in the bulk of the paper does not lead to a well defined v. 
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Figure 17: The L2 norm of the constraint for various values of the resolution as a function of 
the field theory time. 
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